Modeling Lithospheric Deformation Using a Compressible Visco?Elasto?Viscoplastic Rheology and the Effective Viscosity Approach

Deformations of the colder regions lithosphere mainly occur in frictional regime. In geodynamic models, plastic deformations are often highly localized (shear bands) and used as proxies for faults. However, capturing generation evolution shear bands models is troublesome. Indeed, mesh dependency lack convergence affect, to some extent, results models. Here we extend most common plasticity implementation codes (effective viscosity approach [EVA]) include combined effects elastic compressibility, dilatancy, strain softening, viscoplasticity. The latter acts a regularization that cures known issues related plasticity. Using regularized based on M2Di MATLAB routines, show volumetric elasto-plastic can significantly impact crustal-scale banding. We also artificial overstress caused by viscoplasticity be mitigated employing power-law Furthermore, demonstrate algorithms geodynamics (based EVA) accurate those obtained with typically engineering (return mapping consistent tangent operator). Finally, examples long-term tectonic using state-of-the art code MDoodz. They indicate viscoplastic efficiently obtain reliable simulations geodynamics.